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! Abstract 

In this paper, we study the 3D strip packing problem in which we are given a list of 3- 
dimensional boxes and required to pack all of them into a 3-dimensional strip with length 1 
i— i. and width 1 and unlimited height to minimize the height used. Our results are below: i) we 

give an approximation algorithm with asymptotic worst-case ratio 1.69103, which improves 
£^ • the previous best bound of 2 + e by Jansen and Solis-Oba of SODA 2006; ii) we also present 

^ ' an asymptotic PTAS for the case in which all items have square bases. 

! 1 Introduction 

> : 

For packing 2D items into bins or a strip, it is a natural idea to exploit techniques for packing 
lower dimensional (i.e., ID) items. The two-stage packing is particularly well-known: As shown in 
t**"* ' Fig. ^ (a), a bin (or a strip) is divided into shelves and each shelf contains a single layer of items. 
After packing items into shelves, the problem of packing shelves into bins (or a strip) obviously 
becomes the ID bin (or strip) packing problem. The idea originally comes from cutting a large 
board into smaller items efficiently JU]; one can cut the board only in two stages, i.e., cutting 
O ■ horizontally first and then vertically. 

It should be noted that many existing 2D packing algorithms jSJ 0] are based on this two- 
stage packing. In 2002, Caprara [3] established the relation between 2D Bin Packing (2BP) and 
2D Shelf Bin Packing (2SBP). Namely the maximum ratio between the optimal cost for 2SBP 
and that for 2BP is equal to = 1.691... which is the well-known approximation factor of 
the Harmonic algorithm for ID Bin Packing 19 . (A similar relation between 2D Strip Packing 
(2SP) and 2D Shelf Strip Packing (2SSP) was also established by Csirik and Woeginger 8 .) As an 
important byproduct, Caprara also showed an approximation algorithm for 2BP whose asymptotic 
worst-case ratio is arbitrarily close to T^, which first broke the barrier of two for the upper bound 
on the approximability of this problem. 

Our contribution This paper extends the two-stage packing into the 3D Strip Packing (3SP) 
and obtains an approximation algorithm whose asymptotic worst-case ratio is arbitrarily close to 
Too. Our model is standard (see Section |2] for details) and the previous best bound is 2 + e by 
Jansen and Solis-Oba [Tl|- We also show that there is an APTAS for the special case in which all 
items have square bases. 
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Our algorithms use a segment as shown in Fig. ^ (b) instead of a shelf in the 2D case. For 
packing items (whose three sides are all at most 1.0) into a segment, we first divide a segment 
into slips and pack the items into slips by the next-fit (NF) algorithm. The key idea is to make the 
height c of each segment sufficiently large (within a constant), which effectively kills the inefficiency 
of the algorithm for the vertical direction in the sense that the unused space at the top of the 
segment is relatively small. After packing items into segments of the fixed height (=c) and fixed 
length (=1.0), we can obviously use a one-dimensional bin packing algorithm to pack segments. 
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Figure 1: Shelves and segments 



Previous results: On 3D Strip Packing, Li and Cheng presented the first approximation 
algorithm with asymptotic worst case ratio 3.25. Two years later, they gave an online algorithm 
for the problem with asymptotic worst-case (competitive) ratio arbitrarily close to (1. 69103. ..) 2 ~ 
2.8596 ^B]- Then Miyazawa and Wakabayashi |21| 122j improved the asymptotic worst-case ratio 
to 2.67 and 2.64. Very recently, Jansen and Solis-Oba [Tl] improved the asymptotic worst-case 
ratio to 2 + e. 

On 2D Strip Packing, Coffman et al. jS] presented algorithms based on NFDH (Next Fit 
Decreasing Height) and FFDH (First Fit Decreasing Height), and showed that the respective 
asymptotic worst-case ratios are 2 and 1.7. Golan |13j and Baker et al. improved the bound to 
4/3 and 5/4, respectively. Using linear programming and randomization techniques, an asymptotic 
fully polynomial time approximation schemes (AFPTAS) was given by Kenyon and Remila |16j . 

On 2D Bin Packing, in 1982, Chung, Garey and Johnson presented an approximation 
algorithm with asymptotic worst-case ratio at most 2.125. Caprara [3| improved the upper bound 
to 1.6910.... On the other hand, Bansal et al. |2j showed that the 2D bin packing problem does 
not admit an APTAS. Chlebik and Chlebikova [4] further gave an explicit lower bound 1 + 2196 • 
Since the 2D bin packing problem is a special case of the 3D strip packing problem, the lower 
bound holds for 3D strip packing too. 

2 Problems and Notations 

Our model is exactly the same as [T3|- Given an input list L of n three-dimensional boxes, in 
which each box has length, width and height at most 1 respectively, 3SP is to pack all boxes 
into a 3D strip (rectangular parallelepiped) of width 1, length 1 and minimum height, so that the 
boxes do not overlap. In this paper we consider the orthogonal version of the problem without 
rotations, i.e., the boxes must be packed so that their faces are parallel to the faces of the strip 
and the boxes are oriented and cannot be rotated. The problem is obviously NP-hard. For 
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approximation algorithms, we use the standard measure to evaluate them, i.e., the worst-case 
ratio. In this paper, we consider the asymptotic worst-case ratio. Given an input list L and an 
approximation algorithm A, we denote by OPT(L) and A(L), respectively, the height used by an 
optimal algorithm and the height by algorithm A for list L. The asymptotic worst-case ratio 
of algorithm A is defined by 

R% = lim sup{A(L)/OPT(L)\OPT(L) =n}. 

n—*oo £ 

3 Basic tools for algorithms and their analysis 

Fractional Bin Packing (FBP). The continuous version of bin packing plays an important role 
in designing an asymptotic PTAS 2(2 E3- We first give its definition and some properties. Given 
an instance I of one dimensional bin packing, suppose that there are p distinct sizes of the items 
in /, where p is a constant. Let si > S2 > ■■■ > s p be the distinct item sizes and rij be the 
number of items of size Sj for j = 1, ... ,p. A feasible pattern is a vector v = (v\, ... ,v p ) such 
that Y2j=i v j s j — 1) i- e -> an items in a feasible pattern would fit in one bin. Let v denote the 
collection of all feasible patterns for / and v % = (v\, . . . ,v p ) denote the i-th pattern in v, where 
Vj is the number of items of size Sj in the i-th pattern. We further denote X{ to be the number 
of bins being needed for packing the i-th feasible pattern in v. If we allow xi to be a fractional 
number, then the problem becomes the fractional bin packing problem (FBP) and corresponds to 
the following Linear Program (LP): 



Min Erfei/ 3 * 

Xi > 0, v % & v. 



The LP dual of Q is given as follows: 

max. Y?j=i nj^j 
s.t. Ei=i^i<l> v l €v (2) 

VTj > 0, j = l,...,p. 

Optimal values for Q and (J2J) coincide and the following important lemma is due to [H], 

Lemma 1 There exists an optimal solution ir* of {^) such that ir* > ix\ > . . . > tt* (recalling 
s\ > s 2 > ■ ■ ■ > s p ). 

The following lemma [HIES]) says that the optimal values for BP and FBP are almost equal. 

Lemma 2 For any bin packing instance I and for any e > 0, we have OPTbp(I) < (1 + 
e)OPT FB p{I) + 0{e~ 2 ), where OPT FBP {I) is the optimal value for FBP. 

Harmonic algorithm. The Harmonic algorithm was introduced by Lee and Lee ^S]- Given a 
(one-dimensional) bin packing instance I and an integer k > 0, we say an item i belongs to type t 
if its size Sj € (fTT> j] f° r * = 1> • • • ; ^ — 1 an d to type k if Sj € (0, where A: is a constant. Then 
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the Harmonic algorithm packs items of different types into different bins. During packing, if the 
current item of type t does not fit in the corresponding bin, then the algorithm closes the bin and 
opens a new one. Given an item of size x, we define a weighting function fk(x) as follows: 

r i, if ^ < x < i with 1 < t < k, 

JkK > ~\ if < x < ±. 

Let t\ = 1, = ti(ti + 1) for i > 1. For a positive integer k, let m(fc) be the integer such that 
t m (k) < k < t m ( fc )+i. T k = Ya^i % + t m( fc )+ i ' ife=T- Note tnat To ° = lim fc^oo7fc « 1.69103 
The weighting function /^(x) satisfies the following property (see |19| k 

Lemma 3 For each sequence x±, . . . , x m with Xi E (0, 1] and Y%=i x i — 1» 

i=l 



NFDH packing. NFDH was first proposed by Meir and Moser [20] for packing a set of squares 
into a rectangular bin, but NFDH packing can also be applied to pack rectangles. It simply 
works as follows. First sort all rectangles in non-increasing order of their heights. Then pack 
them into the bin level by level and in each level we use the Next Fit (NF) algorithm. If a level 
cannot accommodate the current rectangle, then we close it (will never be used again) and open 
a new one. (see Figure 121(c)). Note that NFDH packing can be extended for multidimensional 
packing |2ffll2]. 

4 New upper bounds for 3D strip packing 

We call our algorithm 3D Segment Strip Packing (3SSP). 
4.1 Algorithm 3SSP 

Given an item R = (l,w,h), where I, w and h are its length, width and height respectively, we 
may use l(R), w(R) and h(R) to denote the three parameters as well. Algorithm 3SSP has the 
following 3 main steps. 

1 . Divide all items into k groups G\ , G2 , ■ ■ ■ , Gk such that those in G{ have their lengths in 
range (^rr, j], where k is a constant. 

2. Sort all Gj-items by their width such that Gi = (R\,R2, ■ ■ ■ ,R ni ) and w(Ri) > 10(^2) > 
• • • > w(R rH ), where n« is the number of items in group Gi for 1 < i < k. Then pack all 
items in Gi, for 1 < i < k, into segments by algorithm GNF (given later). For i = k, pack 
all items in Gk into segments by algorithm GNFDH (given later). 

3. When all items in group Gi have been packed into segments, for 1 < % < k, then regard all 
segments as one dimensional items (Only their width is considered) and call an asymptotic 
PTAS for one-dimensional bin packing (e.g. |15| I26j) to pack these segments. 



4 



In the following we give the procedures to pack 3D items into segments, which are the cores 
of algorithm 3SSP. We deal with G% items (1 < % < k) and G k items separately. 

GNF: Consider d items (1 < % < k). Given G { = R 2 , . . . , R n J such that w(R l ) > w(R 2 ) > 
■■■>w(R ni ). 

1. Open a new segment of size (1, w y , c), where w y <— w(R\). 

2. Divide this segment into i pieces of slips of sizes (j,w y ,c), as shown in Figure H (a), then 
without considering their widths and lengths, pack items into these slips by Next Fit. (see 
Figure [1(b)). 

3. If there are remaining items, re-index them and go to Step 1. 
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Figure 2: GNF for G3; the projection of GNFDH packing without considering widths 



GNFDH: Given G k = (R 1 ,R 2 ,..., R Tlk ) such that w{R x ) > w{R 2 ) >■■■> w{R rik ). 

1. Open a new segment with size (1, w y , c), where w y <— w{R\). 

2. Find a maximal index j such that R\, R 2 , . . . , Rj can be placed into the segment by NFDH 
without considering their widths. Pack the j items by NFDH. (See Figure IU(c)) 

3. Re-index the remaining items in Gk (if any), go to Step 1. 
4.2 Analysis of the algorithm 

In the algorithm analysis, dual feasible functions by Fekete and Scheper [S] play a crucial role. 
(Similar notions are used as weighting functions ^2 El El I2H ES]) Suppose that a function 
/ : [0, 1] — > [0, 1] satisfies YllLi f( x i) — 1 f° r an Y sequence x\, . . . , x m such that YliLi x i — 1 an< ^ 
X{ € [0, 1]. Then / is called a dual feasible function. Here are two specific examples: Let it = 
(tti, ...,vf p ) be a feasible solution of ^ (dual LP for FBP in Section^ satisfying the requirement 
of Lemma^and let 7f p +i := 0, sq := 1 and s p +i := 0. Define a new function g by 

g(0) = 0, and g(x) = ffj, for x G [sj,Sj_i). 

The other example is /& defined in Section |H1 
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Lemma 4 Both g(x) and are dual feasible functions. 

Using these two functions, we define the modified volume W(R) of an item R = (I, w, h) as 

W(R) = f k (l) ■ g{w) ■ h. 

The total modified volume of the input list L of items is W(L) = J2neL W(R). 

We need one more lemma regarding dual feasible functions and 2D packing: let (l\, w^i), • • • ; (fmi ^m) 
be 2D items which can be packed into a square bin of size (1,1), and f\ and fi be dual feasible 
functions. Then we have the following lemma (see for the proof), which is important for 
bounding the total modified volume. 

Lemma 5 YZi h(k)h(™i) < 

Now, we are ready to prove the upper bound for the approximability of our algorithm 3SSP. Let 
I(L) be the 1-dimensional item list obtained after Step 2 of 3SSP, i.e., the list of the widths of 
the segments. Recall that c is the height of the segment and k is the parameter of the Harmonic 
algorithm. Let OPTbp{I{L)) be the optimal cost of 1-dimensional bin packing for the list I(L) 
and OPT(L) be the optimal cost for 3D Strip Packing for the list L. Our goal is thus to prove 
the following theorem. 

Theorem 1 For anye>0,c- OPT BP (I(L)) < ^-(1 + e)T k OPT(L) + 0(cke~ 2 ). 

Since we employ some APTAS for packing I(L), algorithm 3SSP achieves the cost arbitrarily close 
to c-OPTbp(I(L)) in the asymptotic case. It shows that the asymptotic worst-case ratio of 3SSP 
is at most -zr(l + e )^fc f° r anv given e > 0, which tends to as e — > and the constants c and 
k take sufficiently large integers. 

The basic idea of the proof is to establish the relation of the left and right-hand sides of the 
inequalities in the theorem to the total modified volume. Recall that 3SSP uses different segments 
for each G{. A segment is called type i if it contains G{ items. For q = 1, . . . , k, let m g be the 
number of segments of type q and w\ the width of the z-th segment of type q , where 1 < i < m q . 
By algorithm 3SSP, we have 

w\>wl>--->wl q . (3) 
Noting that g(-) is the function defined in Subsection 14.21 for instance I(L), by LemmaG] we have 

g(wl)>g{wl)>--->g(wl iq ). (4) 

Let G q := Y^=i 9( w i) denote the total modified width of the segments of type q. Now, we give a 
lower bound for W(L). For convenience, we define = for all g's. 

Lemma 6 The total modified volume W(L) > (c — 1) X^=i G q ~ 

Proof. Let Sf be the i-th segment of type q and L q be the set of all boxes in Sf. We first prove 
that 

W(Lj)>(c-l)g(wl +1 ). (5) 
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Case 1. q / k. If i = m q , W{L\) > (c - l)g(w q l+1 ) since sKU+i) = and ^(4) > 0. 
Otherwise, i.e., 1 < i < m q , by GNF packing, we have the height packed in every slip of Sf is at 
least (c — 1) and every box in Sf has length in , |] and width at least wf +1 (i.e., the width 
of the next segment). Remember fk{x) = | where x in (^j, |]. By the total weight of the 
boxes in every slip is at least 

1 c — 1 

- • g{w q +l ) • (c - l) = g{w q i+l ). 

Since there are q slips in segment Sf, 

W(Lf) >Q~- = (c - 

Case 2. q = k. If i = m q , W(Lf) > (c — l)g(wf, 1 ) still holds. Otherwise, i.e., 1 < i < m q . 
Consider GNFDH packing for items of type k. Assume there are I levels in Sf, and their heights 
are hi, h%, . . . , hi, respectively. Set hi + i = c — Y^,j=i hi- By NFDH packing, if i < m q , we have 
hi > h% > . . . > h[ + i and every box in has width at least length at most jk Hence the 

total sum of lengths in every level is at least 1 — \ . Remember fk(x) = where x in [0, V\. And 
every box in the j-th. level of Sf has height at least hj+i. So the total weight in the j-th level is 
at least 

k 

hj+i ■ g(wf +1 ) -]—[ l ( R ) ^ %+i • 9(wf+i). 
R 

Since Y2j=i hj > c — 1, the total weight of L\ is 

I 

W{L\) > > (c - 

So the inequality © holds. 

Since there are k types of segments and in every type q there are m q segments, so 

k m q 
q=l i=l 

By ©, 

W(L) > (c-l)ELiE£i^«i) 

= (c-i)Eh(G'-sK)) 

> (c-l)£j =1 G*-cfc. 
The last inequality follows directly from g(w q ) < 1. □ 

Next we give an upper bound for total modified volume W(L). 
Lemma 7 Given any input list L over [0, l] 3 , the total modified volume W(L) < T k OPT(L). 
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Figure 3: Cutting an optimal packing and a layer projection on xy-plane 



Proof. Consider an optimal packing for an input list L. For each item of L we draw two horizontal 
planes at its bottom and top, shown as Figure These planes cut the optimal packing into layers 
such that all items (may be part of the original items) in a layer have the same height. Then we 
can see that each layer is associated with a feasible packing on a square bin of (1, 1) by ignoring the 
heights. Assume that after cutting, totally, there are I layers and their heights are 81,82, ■■■ ,8i, 
respectively. By Lemmas 01 and HJ we have X^eS fk( x ) < Tj. and XxeS g{x) < 1 for any list S 
with YlxeS x — Since in the i-th layer, every item has height 8i, by Lemma [SJ the total weight 
of all items in the i-th layer is at most 

T k x 1 x Si. 

Since 

I 

^8 l = OPT(L), 
i=i 

then 

I 

W(L)<Y J SiT k = T k -OPT(L). 

i=l 

□ 

Now it is straightforward to prove Theorem^ 
(Proof of Theorem ^) By Lemma |5J we have 

OPT BP {I(L)) < (1 + e)OPT FBP (I(L)) + 0( £ - 2 ). 
By the duality of FBP and the dual FBP, as used in |2] , we have 

k 

OPT FBP {I(L))=Y,G q , 

q=l 

where G q := Y2iLi 9( w i) denotes the overall modified width of the segments of type q. By Lemmas 
©andEl ^ 



c-OPT BP (I(L)) < c(l + e)OPT FB p(I(L))+0{ce- 2 ) 



= c(l + e)Y,G q + 0(ce- 2 ) 

.W{L)+ck 2 
< c(l + e) — ^— ■ y0(c(7 2 



< 



c- 1 



c- 1 

(1 + e)T k OPT(L) + 0(cfce- 2 ) 



□ 



Remark. Our algorithm can also be applied to the parametric case in which the boxes have 
bounded length (or width). Then by Theorem ^ the asymptotic worst-case ratio in the parametric 
case that all boxes have width or length bounded from above by a is stated in the following table, 



which is better than the previous parametric ratio RT ara in [2 



a 6 


(ii] 




f 1 -1 


( 1 X l 

v 5' 4J 


poo 

n 3SSP 


1.691... 


1.423... 


1.302... 


1.234... 


^para 


3.050... 


2.028... 


1.684... 


1.511... 



5 APTAS for packing items with square bases 

In this section, by combining the techniques for 2D strip packing ^B] and 2D bin packing [2], we 
give an APTAS for the case that the boxes have square bases (bottoms). 
The standard ideas in our scheme are below: 

• Create a gap between large items and small items such that the items fall into the gap do 
not affect the packing significantly. 

• Pack large items in the way similar to 2D strip packing JH] and pack the other items by 
NFDH [13 HE]. 

We use a multidimensional version of NFDH in |2()1 12"]. called MNFDH, to pack items with small 
base sizes into a 3D bin or a strip. The lemma below can be obtained directly from |20| (see also 

0HH). 

Lemma 8 Let I be a set of 3D boxes with base sides at most 5 and height at most 1. Consider 
the MNFDH packing applied to I. If MNFDH cannot place more boxes from I into a bin of size 
(a,b,c), then either all boxes of I has been packed into the bin or the total packed volume in the 
bin is at least (a — 5)(b — S)(c — 1). 

Given any feasible 3D strip packing of height h, we can scan a plane parallel to the square base of 
the strip from the bottom to the top to obtain a vector x = {x\, . . . , x q ) such that Y2t=i x i = ^> 
where q is the number of patterns to pack all squares induced from the input list into a unit square 
bin and X{ is the height of pattern i. 

Definition of S(K,5). If an input set I has a constant number of different sizes, say K, and all 
the base sides are at least 5, where 5 is a constant, then we define this problem as Restricted 3D 
strip packing with square bases, denoted by S(K,8). 
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Lemma 9 J$ The number of all feasible patterns of packing the square items, induced from an 
instance of S(K,5), into a unit square bin is a constant. 

Lemma 10 S(K , 5) can be solved within OPT + K in polynomial time of n, where OPT is the 
optimal cost for S (K, 8) and n is the input size. 

The proof is put to the appendix. 

Lemma 11 Assume the input set I contains boxes with base sides at least 5. Then for any K > 0, 
we can get a solution within OPT{I)(l + -gr^) + K in polynomial time for packing I into the strip. 

The proof is put to the appendix. 

Asymptotic PTAS Using the similar techniques as in we present an APTAS. Given an 
input set / and any e > 0, our packing is as follows. 

1. Let Wj be the base side length of item j. Define Mj = {j : Wj € [e 2l+1_1 , e 2 * -1 )} for 
i = 1, r + 1, where r = |~l/e] . 

2. Set M := Mj for some index 1 < i < r satisfying Vol(Mi) < eVol(I) (such a set M, must 
exist), where Vol(X) is the total volume of items in X. Define the set of large items as 
L = {j : Wj > e 2 ' -1 } and the set of small items as S = {j : Wj < e 21+1_1 }. 

3. Set K = |T/(e5 2 )] and round all items in L up into K distinct sizes, 5 = e 2 * -1 . Then call 
the algorithm in Lemma ^2 to get an almost optimal solution. 

4. Partition the unused space in the current strip into cuboid regions and use MNFDH to pack 
as many squares in S as possible into the free space. Let S C S denote the subset of the 
remaining small items that could not be packed {S could possibly be empty). 

5. Use MNFDH to pack M U S at the top of the current packing in the strip. 

Theorem 2 Given an input set I of 3D boxes with square bases, A(I) < (l + 12e)OPT(I) + 0(K), 
where A(I) is the height used by our algorithm and K = e~°( 2 ' ). 

Proof. (Sketch.) Our argument is similar as After Step 4, there are two cases. 
Case 1. S' is not empty. Then by the proof in Section 3.4 of [2J, 

A(I) < Vol(I)/(l - 6e) + O(K) < (1 + \2e)OPT{I) + O(K). 
The last step follows by assuming without loss of generality that e < 1/12. 

Case 2. S' is empty. Set K = l/(e5 2 ) = e' ^ \ where 5 = e 2 *" 1 in Step 2. By Lemma ITT1 

A(L US)= A(L) < (1 + -^)OPT(I) + K < (1 + e)OPT(I) + 0{e-°^^). (6) 

Next, we consider the cost of packing M by MNFDH. Since the base size of each item in M is at 
most e, by Lemma |H1 

A(M) < Vol(M)/(l - 2e) + 1 < eOPT(I)/(l - 2e) + 1. (7) 

Combining ® and ©, A{I) < (1 + 3e)OPT(I) + 0(e-°^~^). 

Finally we want to note that each step in our algorithm takes polynomial time of n since e is 
a constant. □ 
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6 Conclusions 



In this paper, we present new asymptotic upper bounds for the 3D strip packing problems. Our 
results give a possible way to apply the approaches for 1- and 2-dimensional bin packing to 3- 
dimensional strip packing. It might be interesting to see if the idea can be used to tackle higher 
dimensional strip packing in the general case. Regrading the special case that items have square 
bases, with the technique in the previous work on 2D bin packing and 2D strip packing an APTAS 
is easily achieved. Such an approach can also be extended to multidimensional strip packing. 
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Appendix 



Proof of Lemma 1101 

Proof. Our idea is similar with the one in [T2|. First we consider the following LP, where q is the 
number of all feasible patterns of packing the squares induced from an instance of S(K, 5) into 
a unit square bin, a^- is the number of type j items in pattern i and (3j is the sum of heights of 
type j items for 1 < j < K, x% is the height of pattern i. 

Min. EUiXi 

Xi > 0, 1 < « < q. 

By Lemma El g is a constant related to K and 5. So, the above LP can be solved in polynomial 
time of n, where n is the input instance size of S(K, 5). Let x* = (x*, x^) be an optimal solution 
of the above LP. By some linear programming property, there are at most K non-zero x*'s. Up to 
renaming, we assume the non-zero coordinates are x"[, ...,x* K . We construct a packing of S(K,5) 
in the following way. 

We fill in the strip bottom-up, taking each pattern in turn. Let x* > be the current pattern. 
Pattern j will be used between level lj = (xj + 1) + • • • + (x*_ 1 + 1) and level lj + \ = lj + x* + 1 
(initially l\ = 0). For each i such that / 0, we draw a.ij cuboids of base size Wi going from 
level lj to level where Wi is the base side length of type i item. After this is done for all j's, 
we take all the cuboids of width W{ one by one in some arbitrary order, and fill them in with the 
boxes of base size Wi in a greedy manner (some small amount of space may be wasted on top of 
each column). 

Since every box has its height at most 1, all boxes can be packed in the above way. Moreover 
x* is a lower bound of the optimal value for S(K, 5). Hence we have this lemma. □ 

Proof of Lemma 1111 

Proof. The algorithm has 3 steps: 

Stacking. Sort the n boxes in non-increasing order of base sizes and stack up them one by one 
to get a stack of height H. And define K — 1 threshold boxes, where a box is a threshold if 
its interior or bottom base intersects some plane z = ^jp-, for 1 < i < K — 1. 

Grouping and rounding. The threshold boxes divide the remaining boxes into K groups. The 
base sizes of the boxes in the first group are rounded up to 1, and the base sizes of the 
boxes in each subsequent group are rounded to the base size of the threshold box below 
their group. This defines oiS(K, 5). 

Packing. Apply the approach of Lemma ITU1 to I sup and output the packing. 

To analyze the algorithm, we use the argument of Kenyon and Remila [TS|. Consider two 
instances I- n j and I' sup derived from the stack built in stacking step. The two instances are 
obtained by first cutting the threshold boxes using the planes z = ijf, then considering the K 
subsequent groups of boxes in turn (where each group now has cumulative height exactly H/K); 
to define I' sup , we round the base sizes in each group up to the largest base size of the group (up 
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to 1 for the first group); to define I' in f, we round the base sizes in each group down to the largest 
base size of the next group (down to for the last group). Let lin(I) be the solution of the above 
linear programming for instance I. It is easy to see that 

lin(l' inf ) < lin(I) < lin(I sup ) < lin(l' sup ). 

Moreover, 

lin(l' sup )<lin(l' inf ) + H/K. 

Since OPT (I) > lin(I) and OPT (I) > 5 2 H, the height used by our packing is at most lin(I sup ) + 
K <lin{I' inf )+H/K + K <OPT{I)+K + H/K < OPT{I){\ + ^) + K. □ 



14 



